Spin polarisabilities of the nucleon at NLO in the chiral expansion 
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Abstract 

We present a calculation of the fourth-order (NLO) contribution to spin- 
dependent Compton scattering in heavy-baryon chiral perturbation theory, 
and we give results for the four spin polarisabilities. No low-energy constants, 
except for the anomalous magnetic moments of the nucleon, enter at this 
order. For forward scattering the fourth-order piece of the spin polar isability 
of the proton turns out to be almost twice the size of the leading piece, 
with the opposite sign. This leads to the conclusion that no prediction can 
currently be made for this quantity. For backward scattering the fourth-order 
contribution is much smaller than the third-order piece which is dominated 
by the anomalous scattering, and so cannot explain the discrepancy between 
the CPT result and the current best experimental determination. 
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Compton scattering from the nucleon has recently been the subject of much work, both 
experimental and theoretical. For the case of unpolarised protons the experimental ampli- 
tude is well determined, and in good agreement with the results of heavy-baryon chiral per- 
turbation theory (HBCPT). However the situation with regard to scattering from polarised 
targets is less satisfactory, not least because until very recently no direct measurements of 
polarised Compton scattering had been attempted. 

The usual notation for spin-dependent pieces of the scattering amplitude in the Breit 
frame is, for incoming real photons of energy u) and momentum q to outgoing real photons 
of the same energy energy and momentum q', 

T = e^e" 

= icr ■ (e f x e) A 3 (u, 0) + icr • (q x q) e' ■ e A^oj, 6) 
+ (icr ■ (e f x q) e ■ q' — icr • (e x q) e • q) A 5 (uj, 6) 

+ (icr ■ (e' x q') e • q' — icr ■ (e x q) e' • q) A e (uj, 9) 

+icr • (q x q) e' ■ q e ■ q' Af(u, 6) + spin independent pieces (1) 

where hats indicate unit vectors. By crossing symmetry the functions A{ are odd in to. The 
leading pieces in an expansion in powers of u are given by low-energy theorems |]J , and the 
next terms contain the spin polarisabilities jf. 

A 3 (uj,6) = -^ r (Q(Q + 2k) - (Q + /t) 2 cos#) + 4vrc«; 3 (7i + 75 cos#) 



A' 



8m* 

2 

At(uj, 9) = - f4(Q + k) 2 + 4vrcA 2 + O(ou 5 ) 
2 

A 5 {u, 6) = ^{Q + k) 2 + 47T^ 3 74 + 0(lu 5 ) 
2 

M"> 0) = + K ) + 47r ^ 3 + °( wS ) 

A 7 (co,9) = 0(cu 5 ) (2) 

where the charge of nucleon is Q = (1 + t$)/2 and its anomalous magnetic moment is 
k = (k s + k v t 3 )/2. Only four of the polarisabilities are independent since three are related 
by 75 + 72 + 274 = 0. The polarisabilities are also isospin dependent. 

None of the polarisabilities has yet been measured directly. The best estimates that exist 
at present are for forward scattering (where only A 3 contributes). The quantity 47r/ 2 (0) is 
defined as dA 3 (u, 0)/do; at uj = 0, and according to the LET it depends only on n 2 [|TJ. 
The relevant polarisability is 70 = 71 + 75, which is related via a dispersion relation to 
measurements at energies above the threshold for pion production, ujq: 

1 p a_(u) -cr+(u) 
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where a± are the parallel and antiparallel cross-sections for photon absorption; the related 
sum rule for the model-independent piece, due to Gerasimov, Drell and Hearn ||, has the 
same form except that 1/uj replaces l/tu 3 . 

Before direct data existed, the relevant cross-sections were estimated from multipole 
analyses of pion electroproduction experiments . These showed significant discrepancies 
between the LET and the GDH sum rule for the difference of /2(0) for the proton and 
neutron, though the sum was in good agreement. Indeed even the sign of the difference was 
different. More recently, measurements have been made with MAMI at Mainz, for photon 
energies between 200 and 800 MeV; the range will be extended downward to 140 MeV, and 
a future experiment at Bonn will extend it upwards to 3 GeV |J. The preliminary data 
from MAMI H suggest a continuing discrepancy between the LET and the sum rule for the 
proton, though a smaller one than given by the multipole analysis. The most recent analysis 
using electroproduction data, which pays particular attention to the threshold region, also 
reduces the discrepancy somewhat |J. 

The MAMI data does not currently go low enough in energy to give a reliable result for 
the spin polarisability, 70 . However electroproduction data have also been used to extract 
this quantity; Sandorfi et al. @]find7o = — l-3x 10~ 4 fm 4 and 7q = — 0.4 xl0~ 4 fm 4 , while the 
more recent analysis of Drechsel et al. |7|] gives a rather smaller value of 70 = —0.6 x 10~ 4 fm 4 . 
(We shall use units of 10 -4 fm 4 for polarisabilities from now on.) 

The spin polarisability has also been calculated in the framework of HBCPT: at lowest 
(third) order in the chiral expansion this gives 70 = a em g\/ (24ii 2 f%ml) = 4.51 for both 
proton and neutron, where the entire contribution comes from nN loops. The effect of the 
A enters in counter-terms at fifth order in standard HBCPT, and has been estimated to 
be so large as to change the sign ||. The calculation has also been done in an extension 
of HBCPT with an explicit A by Hemmert et al. ||. They find that the principal effect 
is from the A pole, which contributes —2.4, with the effect of ttA loops being small, —0.2. 
Clearly the next most important contribution is likely to be the fourth-order ttN piece, and 
this is the result which is presented here. The effects of the A at NLO involve unknown 
parameters; one might hope that the loop pieces at least will be small. 

One other combination of the polarisabilities has been estimated from low-energy data 
for Compton scattering from the proton by Tonnison et al. [|l0j, namely that for backward 
scattering, j n = 71 — 75. This quantity is dominated by the anomalous 7r° exchange graph, 
which vanishes for forward scattering, but at third and fourth order there are also pion loop 
contributions. The experimental value is 7^ = —27.1 with experimental and theoretical 
errors of about 10% each. The HBCPT result of Hemmert et al. is —36.7, of which —43.5 
is the anomalous contribution, 4.6 is the ttN piece and 2.2 comes from including the A ||. 

In the current paper we calculate the fourth-order contributions to all four polarisabilities. 
Results for 70 have been previously presented by Ji and collaborators [|ll|] and by ourselves 

To calculate the spin-dependent scattering amplitude, we work in the gauge A = 0, or in 
the language of HBCPT, v ■ e — 0, where v M is the unit vector which defines the nucleon rest 
frame. Note that in this gauge there is no lowest-order coupling of a photon to a nucleon; 
the coupling comes in only at second order. The Feynman vertex consists of two pieces, one 
proportional to the charge current and one to the magnetic moment: 
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Fig. 1: Diagrams which contribute to spin- dependent Compton scattering in the e • v = 
gauge at LO. The solid dots are vertices from CP^ and the open circle is a vertex from C^ 3 \ 



2M 



(Qe • (pi + pa) + 2(Q + «) [5 ■ e, 5 ■ g]) . (4) 



This and all other vertices are taken from the review of Bernard et al. 

At leading (third) order the fixed terms in the amplitudes A3 to ^6 are reproduced, with 
k replaced by its bare value, by the combination of the Born terms and the seagull diagram, 
which has a fixed coefficient in the third-order Lagrangian || . The loop diagrams of Fig. 1 
have contributions of order uj which cancel and so do not affect the LET. However they do 
give contributions to the polarisabilities. 

At NLO, the diagrams which contribute are given in Fig. 2. There can be no seagulls 
at this order; since the NLO pieces of the Ai(u) are of fourth chiral order and are odd in 
u, they will have expansions of the form e 2 u (am^ + bu 2 /m n + . . .). These non-analytic 
powers of m 2 cannot be present in the basic couplings in the Lagrangian, but can only be 
generated from loops. It follows that there are no undetermined low-energy constants in the 
final amplitude. 

The insertion on the nucleon propagator of Figs. 2a, b, h and i needs some explanation. 
Denoting the external nucleon residual momentum by p, the energy v ■ p starts at second 
chiral order with the mass shift and kinetic energy. In contrast the space components of 
p and all components of the loop momentum I and the photon momentum are first order. 
The propagator with an insertion consists of both the second term in the expansion of the 
lowest-order propagator, i/(v ■ I + v ■ p), in powers of v ■ p/v ■ I, and also the insertions from 
£( 2 \ The second-order mass shift and external kinetic energy cancel between the two. 

To calculate the polarisabilities, it is sufficient to work in the Breit frame in which the sum 
of the incoming and outgoing nucleon three-momenta is zero. In the centre-of-mass frame, 
this is not the case, and at fourth chiral order new terms appear in the amplitudes which 
were absent from the general form of the amplitude given by Hemmert et al. ||. They 
can, however, easily be distinguished from the terms in which we are interested because 
they appear to have the wrong crossing symmetry, contributing at order uj 2 and uj a in the 
expansions of the A^. These terms in fact arise from boost corrections to the third-order 
contributions, the uj 2 terms from the LET pieces, and the a; 4 terms from the third-order 
polarisabilities. There is also one other spin-dependent term — we need to remember that 
a under a Lorentz boost, Wigner rotations can generate spin-dependent amplitudes from 
spin-independent ones. This effect is of second order in the chiral expansion, in contrast 
to the other terms in a boost which are of first chiral order. Since the Thomson term for 
spin-independent Compton scattering is of second order in HBCPT, through the Wigner 
rotation it gives rise to a fourth-order contribution to A4 in any but the Breit frame. 

Returning to the Breit frame, only diagrams 2a-h contribute. When the amplitudes are 
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Fig. 2: Diagrams which contribute to spin-dependent forward Compton scattering in the 
e • v — gauge at NLO. The solid dots are vertices from £( 2 \ 



Taylor expanded, there are contributions at order u and uj 3 . The former do not violate 
the LETs, however; they simply provide the pieces proportional to 5k v = —g 2 A m 7 ,M N /A-nfl 
necessary to satisfy the LET's to fourth order. The contributions of the various diagrams 
from Fig. 2 to the polarisabilities are given in Table 1. It can be seen that the requirement 
75 + 72 + 274 = is satisfied, which provides a non-trivial check on the results. 
The loop contributions to the polarisabilities to NLO are then 



7i 

72 
73 
74 

7o 



Qem9 a 

967r 2 / 2 m 2 

C^emQ a 
O^emQ a 



1 - 
1 - 
1 - 



Trm n 

~4M~ N 
Tim 



+ 5r 3 ) 

+ k v + 3(1 + k s )t 3 ) 
(6 + T3 



71 7TI 

-1 + ^(15 + 4/^+4(1 + /^ 

TTTfl 

1 -^rr( 15 + 3 ^ + ( 6 + ^) r 3' 

oM N 



(5) 



Although the subleading pieces have a factor of m n /M N compared with the leading piece, 
the numerical coefficients are often large. The anomalous magnetic moments are k s = —0.12 
and k v = 3.71; with these values the numerical results for the polarisabilities to fourth order 
are 

71 = [-21.3] +4.5 -(2.1 + 1.3 r 3 ) 

72 = 2.3-(3.1 + 0.7r 3 ) 

73 = [10.7] + 1.1 - (0.8 + 0.1r 3 ) 
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Diagram 
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Table 1: Contributions to the polarisabilities from the diagrams of Fig. 2, in units of 
a em g\/(1927im N M n f 2 ). The polarisabilities are defined in Eq. 0. 



7 4 = [-10.7] - 1.1 + (3.9 + 0.5 r 3 ) 

7o = 4.5- (6.9 + 1.5 r 3 ) 

7. = [-42.7] +4.5 + (2.7 -1.1 r 3 ) 
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The term in square brackets, where it exists, is the third-order anomalous contribution. 
(The fourth-order anomalous contribution is proportional to the difference in the incoming 
and outgoing photon energies, which vanishes in the Breit frame.) 

The NLO contributions are disappointingly large, and call the convergence of the expan- 
sion into question. While the fifth-order terms have also been estimated to be large || , this 
is due to physics beyond pion-nucleon loops, namely the contribution of the A. Our results 
show that even in the absence of the A, convergence of HBCPT for the polarisabilities has 
not yet been reached. However the convergence is better for 7^, and so our results do not 
shed light on the large dicrepancy between the experimental and theoretical determinations. 

After we completed this paper the Jiilich group released their results for the fourth-order 
polarisabilities fl4| . They choose not to include the one-particle-reducible diagram of Fig. 2g 
in their definition of the polarisabilities. With the addition of that contribution, their results 
are in agreement with ours. 

JMcG and MCB acknowledge the support of the UK EPSRC. VK held a Commonwealth 
Fellowship while in Manchester. 



APPENDIX A: FULL AMPLITUDE 

The full amplitude in the Breit frame for diagrams 2a-2h are as follows. The notation ti 
is used for the tensor structures which multiply the amplitudes A4. 



T 



9 2 e 2 



' 2 2/-, , fl A3Jo(w,m 2 ) , 2N 

m — uj (1 + cost*) ) + 2ujJ (uj,m 



T h =- 



g 2 e 2 y2m 2 



2m N fl 
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(j' 2 {uj,m 2 ) - J2(0,m 2 ))t 3 + (1 + cosfl) 



t 3 - (u 

2\ 



-UJ) 



<9J 2 (^, m 
du 
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-uJ 2 (u,m 2 )t 5 + uj(t 5 — 2(1 — cos0)t^j J dxJ 2 (xuj,m 



— [u — > — u) 



g 2 e 2 



g 2 e 2 
m N f 2 

rp 

6 2m w /a 



t 3 c<j / dxJ 2 (xuj, m 2 )t 3 — (a; — > — c<j) 
; JO 



7 



# 2 e 2 



T h = - 



(1 - r 3 )-(j 2 (^,m 2 ) - J 2 (0,m 2 ))t 3 - (w -»• -w) 
-u(2(ii v - n s T 3 )(t 3 cos0 - U) + (1 - T 3 )t 6 ) ( dx(l - 2x) J 2 (xu,m 2 ) - (. 
(2(//„ + n s T 3 )(t 3 cos9 + U — t 5 ) + (1 + r 3 )i 6 ) ^ dxJ 2 (xu,m 
^ ^ fdy r d x[((7x - l)(t 6 - t 5 ) + 7(1 - x - y) U ) *5<*' ^ - XVt) 



9 2 e 2 

T 9 = ''A J5 U 

Am N fl 



UJ — > — UJ) 



uj — > —u 



m N fl Jo Jo 
+ (2V(x, y, 9) (xt 6 -xt 5 + (l-x- y)t 4 ) 

-(1 - x - y)(9xy - x - y)t 7 ^uj 
8Jq(uj, m 2 — xyt) 



do) 



, dJ 2 (cu, m 2 — xyt) 



-xy(l - x- y)u 4 V(x, y, 6)t 7 
where uj = (1 — x — y)uj, 



duj 



duj 

— (uj — > — uj) 



(Al) 



/ 9\ // 2 2\ t / 2\ (jJTTt \ 

J 6 {uj,m ) = ^j— )J2{oj,m ) — A^J, 



(A2) 



J (a;,m 2 ), J 2 (uj,m 2 ) and A^ have their usual meanings, prime denotes differentiation with 
respect to m 2 , and 



V(x, y, 6) = (2xy — x — y + 1) cos 6* — x(l — x) — y(l — y). 



(A3) 
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